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ON  HIGHER  ORDER  DYNAMICS  IN  LATTICE-BASED  MODELS  USING 

CHAPMAN-ENSKOG  METHOD 
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Abstract.  In  this  paper,  we  investigate  the  existence  of  higher  order  dynamics  in  lattice-based  models. 
We  have  identified  two  conditions  that  determine  whether  a  model  would  allow  some  Burnett-like  equations 
when  the  Chapman-Enskog  expansion  is  used.  These  two  conditions  are  the  number  of  the  conserved 
quantities  as  well  as  the  space  and  time  discretization.  We  shall  demonstrate  these  conditions  by  discussing 
(1)  pure  diffusion  equation,  and  (2)  hydrodynamic  equations.  While  the  fact  that  diffusion  equation  allows 
the  higher  order  dynamics  can  be  shown  easily,  we  will  illustrate  that  care  must  be  taken  when  deriving 
Burnett-like  equations  for  lattice-based  hydrodynamics  models  using  the  Chapman-Enskog  method. 
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1.  Introduction.  Compared  to  traditional  methods  in  computational  fluid  dynamics  (CFD),  the  lattice- 
based  models  are  simple  and  easy  to  implement  on  computers.  The  advantages  and  disadvantages  of  the 
original  lattice  gas  automata  (LGA)  have  been  well  documented  [1-7].  The  lattice  Boltzmann  equation 
(LBE)  was  later  introduced  to  remove  some  of  the  drawbacks  [8-10].  A  further  simplification  to  the  LBE  is 
achieved  using  the  BGK  procedure  (LBGK)  [11-14]. 

In  lattice-based  models,  it  is  well  established  that  the  Navier-Stokes  equation  can  be  deduced  at  low  order 
expansion  of  Chapman-Enskog  expansion  [15].  Many  authors  further  asserted  that  the  Burnett-like  equation 
could  be  obtained  by  performing  higher  order  using  Chapman-Enskog  expansion  [4,6,7].  The  motivation  of 
this  paper  is  to  carry  out  these  higher  order  Chapman-Enskog  expansion  to  investigate  whether  it  is  consistent 
to  do  so.  We  will  first  study  the  lattice-based  model  for  pure  diffusion  model  [16,17];  and  demonstrate  that 
higher  order  dynamics  is  allowed  in  this  case.  We  will  then  point  out  that  the  Burnett-like  equations  could 
be  derived  for  lattice-based  hydrodynamics  models.  Attention  should  be  paid,  however,  when  the  classic 
Chapman-Enskog  expansion  is  applied  because  of  the  non-commutative  feature  of  cross  derivatives  of  two 
time  scales,  these  derivatives  do  not  exist  in  the  continuous  time  and  space  while  do  exist  in  discrete  velocity 
models  [18].  The  number  of  conserved  quantities  is  also  critical  for  the  existence  of  higher  order  equations. 

2.  High  Order  Dynamics:  Pure  Diffusion.  We  now  consider  the  lattice  BGK  models  for  pure 
diffusion  problems  where  the  only  quantity  conserved  during  the  redistribution  is  the  total  mass.  The 
propagation  step  is  the  same  as  lattice  gas  models  while  the  collision  step  is  just  a  redistribution  of  mass  in 
all  possible  directions.  We  start  with  the  following  evolution  equation  [12], 

(2.1)  fi{x  +  cl,t  +  1)  =  +  u>(f-q{x, t)  -  fi(x, t)) 
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where  fi  is  the  average  population  of  particles  with  velocity  c£(i  =  1,2, ...,  B)  which  belongs  to  a  predeter¬ 
mined  finite  set  and  u>  the  relaxation  parameter  which  satisfies  0  <  u  <  2.  The  local  equilibrium  population 
is  chosen  as  [17], 

(2.2)  f?q(  £,  t)  =  Wip(£,  t),  wt  = 

B  is  the  number  of  particles’  discrete  velocities.  This  is  a  homogeneous  equilibrium  population  in  all  velocity 
directions.  The  macroscopic  density,  denoted  by  p,  is  defined  by: 

B  B 

(2.3)  p(x,  t)-J>  {x,  t)  =  £  f-q  (x,  t) . 

i= 1  *=1 

The  weighting  factor  m*  satisfies  the  normalization  constraint:  =  1.  The  choice  (2.2)  for  the  equilib¬ 

rium  population,  when  used  together  with  (2.1)  and  (2.3),  will  be  shown  to  lead  to  the  diffusion  equation. 
We  consider  models  with  the  particle  velocity  set  in  D  dimension  ( D  =  1,2  and  3).  The  simplest  models 
take  the  velocity  set  of  2D  elements:  D  directions  along  axis  and  D  opposite  directions.  The  rest  particles 
can  also  be  included. 

We  assume  a  weak  deviation  from  the  local  equilibrium  /?9(£,  £), 

(2.4)  fi{x,  t)  =  ftq  (x,t)  +  (x,  t )  +  e2/f )  (x,t)  +  ■■■ 

where  e  is  the  appropriate  Knudsen  number.  The  space  and  time  derivatives  are  expressed  in  terms  of 
multiple-scale  variables  up  to  the  fourth  order  in  time  (see,  for  example,  Huang  [19]), 

(2.5)  da  =  eda 

(2.6)  dt  —  cdtx  +  c*dt2  +  +  e49t4. 

When  the  total  mass  is  conserved,  it  follows  from  (2.1),  (2.2),  (2.3)  and  (2.4)  that, 

B 

(2.7)  £^  =  0,  j  >0. 

i=l 

Using  the  classic  Chapman-Enskog  expansion  and  taking  into  account  of  the  discreteness  of  lattice  model, 
we  obtain  the  first  order  equation  in  e, 

(2.8)  dtlp  =  0. 

The  second  order  equation  is, 

(2.9)  dt2p  -  |p(^  -  1  )daap  =  0. 

The  equations  (2.8)  and  (2.9),  i.e.,  the  dynamical  equations  from  the  two  separated  time  scales  1/e  and  1/e2, 
are  now  reconstituted  to  obtain  the  macro-dynamical  equations  for  the  model.  The  equation  of  diffusion 
equation  is  obtained  from  (2.8)  and  (2.9) 

(2.10)  dtp  =  K2daap 
where  the  diffusivity  k2  is  given  by 

(2.11)  «  = 
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Fig.  2.1.  The  dispersion  relation  (up  to  fourth  order)  versus  k  for  the  DSQ6  model ,  The  open  triangles ,  solid 
triangles,  open  squares ,  solid  squares  and  open  circles  are  numerical  simulations  corresponding  to  w  =  0.75, 1.0, 1.25, 1.5  and 
1.75,  respectively.  The  critical  value  cvcr  is  1.0  for  this  model. 

We  can  also  obtain  higher  order  equations  by  carrying  the  Chapman-Enskog  expansion  further.  We 
derive  the  third  order  equation, 


and  the  fourth  order  equation, 


dt3p  —  0 


dt4p  =  —Aidaappp  —  A2daaaocP • 


The  coefficients  Ai,A2  and  «2  in  (2.10)  for  models  including  rest  particles  are  obtained  after  some  algebraic 
calculations, 

,  v  c?,2 

(2-14) 


(2.15) 


.  c2 .  2  2  1, 

A\  =  o - h  7)^2 

D  a;'*  cj  4 


(2.16) 


The  final  fourth  order  equation  is  the  following  [17], 

(2.17)  &tP  ~  K2 SococP  -A,  & j3 P  A  2 ^aocococP' 

We  note  that  Equation  2.17  is  anisotropic  due  to  the  last  term.  Applying  the  Fourier  transform  exp(-Clt  - 
ikx)  ( k  is  the  wavenumber  and  fl  the  frequency)  to  the  above  equation  in  one-dimensional  space,  we  get  the 
dispersion  relation  which  reads  as, 

(2.18)  —  =  1  +  —  k2, 

K2  «2 


3 


where  ka  =  A\  +  A2  and  k  =  p-. 

Numerical  result  is  given  by  the  Figure  2.1.  The  curves  correspond  to  theoretical  results  k/k2  while  the 
points  correspond  to  numerical  simulations  of  the  lattice  model  presented  above.  Satisfactory  agreements  in 
all  cases  are  achieved.  The  fourth  order  corrections  may  have  effects  in  the  regime  of  large  Knudsen  number, 
i.e.,  large  k  and  small  u).  Equation  2.18  is  valid  only  for  wavevector  along  x  (or  y,z)  axis,  so  is  the  critical 
value  ucr  —  1  for  the  D3Q6  numerical  model  [12]  used  for  Equation  2.1. 

3.  High  Order  Dynamics:  Hydrodynamics.  We  now  turn  our  attention  to  lattice-based  hydrody¬ 
namics  models.  In  the  LGA,  LBE,  and  LBGK  models,  both  the  mass  and  momentum  are  conserved.  The 
common  features  in  these  models  are  discrete  velocity  space  of  particles,  evolution  steps  of  local  interac¬ 
tions  and  neighbor-to-neighbor  propagation  of  moving  particles.  Since  the  principle  of  deriving  large-scale 
equations  is  the  same  and  outlined  in  the  previous  section.  For  the  sake  of  simplicity,  we  use  lattice  BGK 
models  to  illustrate  the  existence  of  high  order  dynamics:  Burnett-like  equations.  In  classic  kinetic  theory, 
Euler,  Navier-Stokes,  Burnett  and  Super-Burnett  equations  constitute  the  successive  approximations  of  the 
Boltzmann  equation  in  the  order  of  Knudsen  number.  Like  in  classic  kinetic  theory,  the  lattice-based  models 
for  hydrodynamics  use  the  Chapman-Enskog  expansion  in  order  to  derive  the  Navier-Stokes  equations.  We 
outline  the  basic  ingredients  of  the  derivation.  The  time  evolution  equation  is  the  same  as  section  2,  except 
that  the  equilibrium  distribution  f^q  contains  not  only  mass,  but  also  momentum, 


(3.1) 


ft9  —  tpP{  1  + 


CiaV'a  (pia^ip  -^'P 


+ 


2cj 


where  cs  is  a  constant.  The  density  p  and  velocity  u  are  defined  by, 


(3.2) 


i=  1 


B 

1 

i= 1 


B 


i= 1 


i=  1 


which  leads  to  the  constraints  on  high  order  corrections 


B 


(3.3)  ^/i  =  0,  =  °»  ■?  >  °- 

i=l  i=l 

The  leading  order  on  e  yields  the  inviscid  fluid  equations, 


(3.4) 

(3.5) 


dtlp  +  da{pua)  =  0 
dtl{pua )  +  dp(puaup)  —  -c2sdap 


and  the  second  order  e2  results  in  the  dissipative  terms, 


(3.6) 

(3.7) 


dt2p  =  0 

dt2(pua)  =  v(dpp{pua)  +  dap{pup)) 


where  v  is  the  shear  viscosity  (v  =  e2(l/u  —  1/2)). 

Now,  in  order  to  obtain  high  order  hydrodynamical  equations  of  the  lattice-based  models,  let  us  look  at 
the  third  order  e3,  the  Taylor  expansion  gives  the  following  equation, 
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dtafiq  +  Ciadafi^  +  d^fj^  +  dt2f. —  ®Mk  + 

—  (dh  t,  +  2ciadtia  +  CiaCipdapS^fi  ^  +  g  (9{1t1t1  +  3Ciadt1t1a  +  3 CiaCipdtl0cj3  +  (kacipCi~fdap-y)f^ 
(3.8)  =  -w/i3). 

Summing  the  underlined  cross  derivative  dtlt 2f°q  in  the  above  equation  over  i,  we  get  a  term, 

9txtM- 

Using  the  first  and  second  order  Equations  3.4-3. 7,  we  obtain  two  different  results, 

(1) .  if  we  first  take  the  derivative  over  f2  then  t\,  we  have, 

dtitiip)  =  0. 

(2) .  Reversely,  we  have, 

dtit2{p)  =  -vda{dpp{pua)  +  dap{pup)) 

It  means  that  the  operators  are  not  commutative, 

9tit2(*)  i1  dt2t!  (•) 


where  •  is  either  p  or  pua. 

Note  that1  the  third  order  macroscopic  equations  can  be  also  obtained  by  the  wavevector  expansion  (see 
for  example,  van  Coervorden  et  al.  [20]).  Even  though  the  above-mentioned  operators  are  not  commutative, 
the  essential  point  in  the  Equation  3.8  is  the  sum  of  the  two  terms.  After  a  tedious  algebraic  calculation,  we 
get  the  third  order  equations, 

c2 

(3.9)  dHp  =  -^dapp(pua) 

(3.10)  dt3(pua )  =  |(^  -  §  +  1  )dapp{p). 

We  check  the  dispersion  relation  up  to  the  third  order  numerically  in  Figure  3.1  (the  curves  are  theoretical 
predictions  with  Equations  3.9-3.10  and  points  numerical  simulations).  Good  agreement  is  obtained. 

Even  higher  order  (fourth  and  up)  dynamics  can  be  obtained  while  tremendous  care  has  to  be  taken 
since  more  non-commutative  operators  are  involved  and  results  will  be  published  elsewhere. 

4.  Concluding  Remarks.  In  this  paper,  we  pointed  out  that  two  conditions  determine  whether  the 
lattice-based  models  could  or  could  not  have  higher  order  dynamics  when  classical  Chapinan-Enskog  ex¬ 
pansion  is  used.  These  conditions  are  number  of  conservation  laws  and  the  space  and  time  discretization. 
The  pure  diffusion  model,  a  system  with  only  one  conserved  quantity,  is  first  presented  to  illustrate  that  the 
higher  order  dynamics  is  allowed.  We  then  turned  our  attention  to  the  lattice-based  hydrodynamics  equa¬ 
tions.  With  more  than  one  conserved  quantities,  we  note  that  special  care  must  be  taken  to  derive  governing 
equations  for  higher  order  dynamics.  After  noting  the  feature  of  no-commutative  cross  time  derivative,  we 
demonstrate  how  Burnett-like  equations  could  be  obtained  for  lattice-based  hydrodynamics  models  using 
the  classic  Chapinan-Enskog  expansion  method.  The  results  reported  in  this  paper  can  be  used  to  analyze 
theoretically  systems  where  hydrodynamic  description  may  break  down,  a  typical  example  is  simulations  of 
the  micro-electronic  mechanical  systems  (MEMS)  [21,22]. 

1Thc  authors  arc  very  grateful  to  the  referee  of  the  Phys.  Rev.  E  for  tills  and  several  other  important  observations. 
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Fig.  3.1.  The  dispersion  relation  (up  to  third  order):  The  speed  of  sound  versus  k  for  the  DIQ5  model,  The  open  triangles, 
solid  triangles ,  open  squares,  solid  squares  and  open  circles  are  numerical  simulations  corresponding  to  u  =  0.75, 1.00,  and 
1.50  while  the  curves  are  theoretical  predictions . 
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